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Abstract 

Neutrino oscillations, biased by the magnetic field, alter the shape of the neutrinosphere in 
a cooling protoneutron star emerging from the supernova collapse. The resulting anisotropy 
in the momentum of outgoing neutrinos can be the origin of the observed proper motions of 
pulsars. The connection between the pulsars velocities and neutrino oscillations results in a 
prediction for the r neutrino mass of m(v T ) ~ 100 eV. 
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It is well known that rotating magnetized neutron stars, pulsars, exhibit rapid proper 
motions 0, |2| characterized by space velocities that range in the hundreds of kilometers per 
second. Born in a supernova explosion, a pulsar may receive a substantial "kick" velocity 
due to the asymmetries in the collapse, explosion, and the neutrino emission affected by 
convection || |J. Evolution of close binary systems may also produce rapidly moving pul- 
sars [Q. Alternatively, it was argued that the pulsar may be accelerated during the first 
few months after the supernova explosion by its electromagnetic radiation, the asymmetry 
resulting from the magnetic dipole moment being inclined to the rotation axis and offset 
from the center of the star. Most of these mechanisms, however, have difficulties explaining 
the magnitudes of pulsar spatial velocities, which can be as great as 1000 km/s and have an 
average value of 450 km/s [|T|. (A recent study J3J shows, however, that a kick velocity of 
this magnitude can be achieved in an asymmetric collapse.) 

In this letter we suggest an explanation for the birth velocities of the neutron stars based 
on asymmetric emission of neutrinos due to neutrino oscillations. We also show that the 
distribution of pulsar velocities can yield information about neutrino masses. 

The basic idea is the following. Neutrinos emitted during the cooling of a protoneutron 
star have total momentum, roughly, 100 times the momentum of the proper motion of 
the pulsar. A 1% anisotropy in the neutrino distribution would result in a "kick" velocity 
consistent with observation. In the dense neutron star an electron neutrino, v e , has a shorter 
mean free path than u^, u T , or any of the antineutrinos. If one of the latter, e. g., z/ r , undergoes 
a resonant oscillation into z/ e , above the r-neutrinosphere but below the e-neutrinosphere, 
it will be absorbed by the medium. Therefore, the effective r-neutrinosphere in this case is 
determined by the point of resonance (Fig [I]). 

It is known that the electromagnetic properties of a neutrino propagating in medium, in a 
longitudinal magnetic field B, are different from those in a vacuum [[5J. The effective neutrino 
self-energy has a contribution proportional to (B ■ k), where k is the neutrino momentum. 
The position of the resonance of the v T — > u e oscillations is affected by the magnetic field and 
depends on the relative orientation of B and k ||. Therefore the effective r-neutrinosphere 
(or, "neutrinosurface" , to be more precise) is not concentric with the electron neutrinosphere 
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Figure 1: If the resonant oscillation v T — > v e takes place between the electron neutrinosphere, 
S e , and that of u^, is T , v T and P e , S^, then the effective r-neutrinosphere coincides with 
the surface of resonance, S^ eS \ The latter is affected by the magnetic field. Therefore, the 
r-neutrinos emitted in different directions come from the regions of different temperatures. 
The resulting anisotropy in the momentum of the outgoing neutrinos can be the origin of 
the pulsar "kick" velocity, whose magnitude depends on the magnetic field. 
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(and, in fact, is not a sphere). The r neutrinos that escape in the direction of magnetic 
field have, therefore, a different temperature from those emitted in the opposite direction. 
They carry away different momenta, thereby creating an asymmetry in good quantitative 
agreement with the data. We will show that this mechanism may be the origin of the proper 
motions of pulsars. 

In a recent analysis of a sample of 99 pulsars, Lyne and Lorimer [IJ concluded that the 
space velocity of pulsars at birth has a mean value of 450 ± 90 km s _1 . Typical pulsars 
have masses in the range from 1.0 M to 1.5 M & , where M & is the solar mass 0. The 
momenta associated with the proper motion of pulsars are therefore of order k p ~ (1.2 ± 
0.4) x 10 41 gcm/s. The energy carried off by neutrinos in a supernova explosion is estimated 
to be ~ 3 x 10 53 erg \W\. Since p u ~ E v (we set h — c — 1), this corresponds to the sum of the 
magnitudes of the neutrino momenta « 10 43 g cm/s. Comparing this with k p , we conclude 
that a few per cent asymmetry in the distribution of the outgoing neutrinos is sufficient to 
give the pulsar a "kick" velocity ~ 450 km s _1 . 

The average energy of the neutrinos emitted during the cooling of the protoneutron star 
depends on the temperature of the neutrinosphere TCfl , defined as the surface where the 
optical depth, dr/\ u , becomes of order 1. Neutrinos of different types have, in general, 
different neutrinospheres because of the difference in the mean free path, \ v oc l/a u . Roughly 
speaking, one can consider the inside of the neutrinosphere to be opaque for neutrinos, 
while the outside is transparent. Electron neutrinos have both charged and neutral current 
interactions; they scatter off electrons, positrons and nucleons, as well as being absorbed |T| 
in the reaction v e n — + e~p. On the other hand, /x and r neutrinos and antineutrinos can only 
scatter elastically, via neutral currents, at MeV temperatures. (The electron antineutrino, z/ e 
can scatter and be absorbed in the process u e p — > e + n, so its opacity starts out close to that 
of v e and then decreases to, roughly, that of and v T as the density of protons decreases.) 

This difference in scattering cross-sections gives rise to roughly an order of magnitude 
difference [[11], [E| in the mean free paths of the v e and the v T in the vicinity of the electron 
neutrinosphere, where the density is 10 11 — 10 12 g/cm 3 . Consequently, the r-neutrinosphere 
lies deeper inside the protoneutron star than the electron neutrinosphere. As a result, the r 
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neutrinos escape from a deeper layer of the star, where the temperature is higher, and their 
mean energy exceeds that of ^ e 's by about 50%. 

We now consider neutrino oscillations. For defmiteness, we discuss two flavors, u e and u T , 
with Am 2 = m 2 (v T ) — m 2 (z/ e ) « m 2 (v T ) ^> m(u e ) « and small mixing. We will also assume 
that u e — > ^ oscillations take place in the Sun, while the transition v T — > z/ e occurs at a 
much higher matter density. We leave a more careful analysis of the neutrino mass matrix 
consistent with all experimental data for future publication |L3] . 

As was shown in Ref. ||, the neutrinos of energy E « k = \k\ propagating in the 
degenerate electron gas of the charge density N e = n e - — n e +, in magnetic field B, undergo 
a resonant oscillation^] if 



where 9 is the neutrino mixing angle in vacuum. 

If the resonant conversion of v T into u e occurs at some point between the r neutrinosphere 
and the electron neutrinosphere, the r neutrinos produced inside the "resonance-sphere" ii. e. 
closer to the center than the point of the resonance), will turn into z/ e 's and will be absorbed 
(thermalized) by the medium. The r-neutrinos produced by the v e — > v T oscillations will 
escape with the energy determined by the temperature at the resonance point. Of course, 
r-neutrinos produced outside the resonance- sphere will also escape, but there is no effective 
mechanism for copious production of r-neutrinos outside neutrinosphere, where the matter 
density is small. The surface of resonance points for the v T v e oscillations becomes the 
effective r neutrinosphere. Since the last term in (0) depends on the relative orientation of k 
and B, the resonant oscillations occur at different densities for the neutrinos going in different 
directions. The corresponding difference in temperatures will result in the asymmetry that 
gives the neutron star a "kick" in the direction of the magnetic field. 

In the absence of the magnetic field, the resonance will occur at a distance r from the 

1 We emphasize the difference between these chirality-preserving oscillations and the spin and flavor 
precession of neutrinos in magnetic field studied, e. g., in Ref. (l4|. The last term in equation (|l|) arises due 
to the effective interaction of neutrinos with the charged particles in the background and does not depend 
on the size of the neutrino magnetic dipole moment. 
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center if 



Am 2 cos26 = 2V2G F N e (r )k. (2) 

If r is to lie between the two neutrinospheres, where the density is ~ 10 11 g cm~ 3 and 
the electron to baryon ratio Y e « 0.1, then, for a small mixing angle, Am 2 has to be of order 
10 4 eV 2 , which corresponds to a r neutrino mass of 100 eV, consistent with the current 
experimental limits ||15|| . This is close to the limit on stable neutrino masses imposed by 
cosmological constraints, m u < 92 eV (fi/ij)) 10' although the mixed dark matter models 



favor a lower mass range [17] for stable neutrinos. Larger values of mass have been considered 
for unstable z/ r 's [|18f| . 

The oscillations are adiabatic as long as the density N e (r Q ) and the magnetic field can 
be considered constant over the oscillation length, l osc , 

( 1 Am 2 \ 1 1 cm . . 

losc~\ — sin 26) » -. 3 

\2tt 2k J sin26 v ; 

For a wide range of mixing angles, l osc is much smaller than the scale on which the density 
changes are noticeable, so the oscillations can, in fact, be treated as adiabatic. 

In the presence of the magnetic field, the condition (|l|) is satisfied at different distances 
from the center, depending on the value of the (k ■ B) term in ([!]). The surface of the 
resonance is, therefore, 

r{4>) — r o + 3 cos (ft, (4) 
where coscf) = (k ■ B)/k and 5 is determined by the equation: 

/s^% (5) 



dr \ 



7T J 



This yields 



, f '3iV e y/ 3 e / dN e (r) e^ e / dN e (r) 
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where /i e pa {?>n 2 N e ) l /' i is the chemical potential of the degenerate (relativistic) electron gas. 

We now estimate the size of the "kick" velocity from the effect just described. As was 
established in the beginning, a few per cent asymmetry in the momentum distribution of 
neutrinos is necessary to explain the observed pulsar velocities. Approximately 1/6 of the 
total energy is carried off by each of the neutrino species UlQ| . Integration over the angles 
in (k ■ B) gives a factor 1/2 if the magnetic field is uniform. (The latter is, of course, a 
simplification.) The asymmetry in the third component of momentum is, therefore, 

Ak = 1 I T\r -5)-T\r + 5) ^ 2 1 dT 

k 6 2 T 4 (r ) ~ 3 T dr ' U 

where we have assumed a black-body radiation luminosity oc T 4 for the effective neutri- 
nosphere. Using the expression (|5|) for 5, we obtain 

Ak e f dT 



k 3tt 2 \'dN e/ 

where rj = fi e /T is the degeneracy parameter of the electrons. Here we used the identity 
(dT I dr) I (dN e / dr) = dT/dN e . To calculate the derivative in @, we use the relation between 
the density and the temperature for the relativistic Fermi gas: 

f d 3 P 1 

e = 2 J (27r)3e(P-")/ r + l- (9) 
Differentiating the right-hand side with respect to T, we obtain dN e /dT w (2/3) T 2 77, for 

large 77. (This asymptotic expression is accurate to three significant digits for 77 > 5.) We 
observe that for a degenerate relativistic electron gas (77 ^> 1) the product T](dT/dN e ) ps 
1.50/T 2 , independent of 77. Finally, the ratio in @ is 

A A; B /3MeV\ 2 / B \ 

— = 0.015 — = 0.01 — - . 10 

k T 2 V T J V3 x 10 14 G/ v ' 

During the cooling stage of the protoneutron star, the r-neutrinos come out with the 
average energy pa 10 MeV pi, which corresponds to the temperature of ~ 3 MeV. We see 
that the observed pulsar velocities, which require Ak/k to be of order 0.01, can be explained 
by the values of B ~ 3 x 10 14 G. 
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The magnetic fields at the surface of the neutron stars are estimated to be of order 
10 12 — 10 13 G H|]. However, a magnetic field inside the pulsar may be as high as 10 16 G |, @. 
The existence of such a strong magnetic field is suggested by the dynamics of formation of 
the neutron stars, by the stability of the poloidal magnetic field outside the pulsar, as well 
as by the fact that the only star whose surface field is well studied, the Sun, has magnetic 
field below the surface which is 10 3 larger than that outside. (This strong magnetic field 
causes the sun spots when it penetrates the surface.) 

Magnetic fields of order 10 16 G inside the neutron star can themselves lead to asymmetric 
neutrino flux by the creation of neutron star analogue of sun spots. This phenomenon 
is disscussed in the context of strongly magnetized neutron starts by R. C. Duncan and 
C. Thompson In such strong magnetic fields, other weak interactions effects also become 
important ||21|| . We emphasize that the magnetic fields relevant for these effects are an order 
of magnitude higher than those considered in the present work. 

It is clear from equation ([TO) that the deformations of the neutrinosphere due to neutrino 
oscillations biased by the magnetic field can result in the asymmetry of the neutrino flux 
necessary to give the pulsar a "kick" velocity consistent with the data. 

In our calculations, we have assumed that the temperature distribution is not affected 
significantly by the position of the resonance. This assumption is well-justified because only 
the neutrinosphere of one out of six (anti-)neutrino species depends on the resonance, while 
the temperature profile is determined by the emission of all six. However, the increase of the 
heat drain in the vicinity of the effective r-neutrinosphere makes the temperature gradient 
larger. This effect further increases the ratio fliTf). We have also neglected the deviations 
from thermal equilibrium in a cooling protoneutron star. A detailed analysis of the heat 
transport equations, though clearly important, is beyond the scope of this letter. 

According to equation fllPf) , the birth velocities of neutron stars are proportional to the 
magnetic field at the effective neutrinosphere. However, the direction of the proper motion 
of a pulsar need not be correlated with its angular momentum, because of the inclination 
of the magnetic field to the rotation axis and its possible offset from the center. This is 
also in good agreement with the observations which find no correlation between the space 
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velocities and the linear polarization of radio emission [22]. In a separate paper |23 



we 



discuss the implications of our mechanism for the correlation between the magnetic fields 
and the magnitudes of pulsar velocities. 

In conclusion, we have described a new explanation for the birth velocities of pulsars, 
which is in good agreement with the observational data. It presupposes, as we have stated, 
a t neutrino mass of about 100 eV for small mixing. (A more complex realization of our 
mechanism may allow for smaller m(z/ r ) |T5| .) The surprising connection between the proper 
motions of pulsars and neutrino oscillations provides a new astrophysical "laboratory" for 
neutrino physics. Measurements of the space velocities of pulsars can yield information about 
the neutrino masses in the range which is currently inaccessible to either solar neutrino, or 
collider experiments. 

We would like to thank A. Dolgov, V. Kaspi, D. Lai, P. Langacker, M. Ruderman, 
R. Sawyer, D. Schramm, D. Seckel, R. Shrock, A. Smirnov and L. Wolfenstein for helpful 
discussions. This work was supported by the U. S. Department of Energy Contract No. 
DE-AC02-76-ERO-3071. 
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